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ABSTRACT 

Context. The prototype of the p Cephei class of pulsating stars, yS Cep, rotates relatively slowly, and yet displays episodic Ha emission. 
Such behaviour is typical of a rapidly rotating, classical Be star. For some time this posed a contradiction to our understanding of 
the Be phenomena as rapid rotation is thought to be a prerequisite for the characteristic emission phases of Be stars. Recent work 
has demonstrated that the Ho- emission is in fact due to a close companion (separation »0.25") of the star This resolves the apparent 
enigma if this close companion is indeed a classical Be star, as has been proposed. 

Aims. We aim to test the hypothesis that this close companion is a valid Be star by determining properties such as its spectral type 
and V sin i. 

Methods. We employed the technique of spectroastrometry to investigate the close binary system. Using the spectroastrometric 
signatures observed, we split the composite binary spectra into its constituent spectra in the B band (4200 - 5000A) and R band 
(6200 - 7000 A). 

Results. The spectroastrometrically split spectra allow us to estimate spectral types of the binary components. We find that the primary 
of the close binary system has a spectral type of B2III and the secondary a spectral type of B5 Ve. From the relationship between mass 
and spectral type, we determine the masses of the binary components to be Mp,.i = 12.6 + 3. 2Mq and M^^^ = 4.4 + 0.7Mq respectively. 
The spectroastrometric data allow some constraint on the orbit, and we suggest a moderate revision to the previously determined orbit. 
We confirm that the primary of the system is a slow rotator (v sin i = 29**1 ^ while the secondary rotates significantly faster, at 
a V sin ; = 230 ± 45 km s"' . 

Conclusions. We show that the close companion to the /3 Cephei primary is certainly a valid classical Be star. It has a spectral type of 
B5Ve and is a relatively fast rotator. We confirm that the /? Cephei system does not contradict our current understanding of classical 
Be stars. 

Key words, binariesxlose - binaries: general - Stars:emission-line,Be - Stars:fundamental parameters - Stars:individual:/? Cephei 



1. Introduction 

The well-studied object /3 Ceph ei is a massive, pulsatin g star 
with a spectral type of early B dHoffleit & Warred [T995h . The 
star is a prototype of a class of early B type giants and sub- 
giants that exhibit rapid radial ve locity, photometric and line pro - 
file variations due to pulsations ( iSterken & Jerzvkiewicz|[T993l) . 
jS Cep lies at a distance of ij=182 +18 pc and has a visual mag- 
nihide of 3.2 (iPerrvman & ESAlll997h . fDonati et al.1 (1200 lb esti- 
mate that the star has a bolometric magnitude of Lhoi - -5.8+0.2 
and a radius of = 6.5 + 1.2Ro. Evolutionary tracks suggest 
that the star's mass is M,-12 + 1 M (T) and that its age is ap- 
proximately 12Myr dPonati et alJl200L) . As such it seems to be 
a fairly typical /3 Cep star, as might be expe cted from the proto- 
type of the class dStankov & Handleiil2005h . 

In fact, jS Cep is a tertiary system. A visual comp anion lies 
at a distance of 13.4"from the primary dHeintzlll978h . The pri- 
mary also has a closer, spectroscopic companion at a separa- 
tion of approximately 0.25"(Gezari et al. 1972). Hereafter, we 
refer to this star as the close companion and the primary of 
the system as /3 Cep. Little is known about the close compan- 
ion. Measurements of the brightness difference between the pri- 
mary and the close companion range from 5.0 magnitudes to 1 .8 
magni tudes (Gezarietal. 1972; Balegaetal. 2002). Catanzaro 
d2008l) indirectly esti mated the vsin; of the clos e companion 
to be «i230 km s ' and lPigulski & Boratynl dl992h estimated its 



mass to be a!8M0. However, usin g the same orbital p arameters 
but with a slightly reduced period Donati et~an d200ll) conclude 
the secondary has a mass of M* - 6.2 + O.3M0. 

Ha emi s sion e manating from /? Cep was first reported 
by iKarpovl dl932h . and has since been found both absent 
and present with a recurring timescale of approximately 
a decade ( Wilson & Se d don 1956; Kaper & Mathias 19951 
iPan'ko & Tarasovl Il997t iNeiner et al.f l200lh . Such episodic 
emission is typical of a classical Be star. However, yS Cep is 
a slow rotator with a rotatio nal period of 12 day s and a v sin ; 
of approximately 25 kms ' (He nrichs et alj|2000l ; iTelting et al.l 
1997). Be stars are typical ly rapid rotators rotating at up to 80 % 
of their break-up velocity dPorter & Riviniusl2003h . As this rapid 
rotation is thought to be integral to the Ha emission of Be stars, 
the Ha err iission of 13 Cep is not consistent with the classical Be 
paradigm. lDonati et al.l (^2001) proposed that the Be behaviour of 
the star is due to a magnetically confined wind leading to shocks 
in the equatorial region (fi Cep does possess an oblique magnetic 
field). However, if this were the case, the Ha emission w ould be 
modulated by the rotation of the star, and it is not ( Schne rret al.l 
2006). Thus the source of the Ha emission has posed a conun- 
drum to our current understanding of Be stars. 

This enigma has been resolved by lSchnerr et alj d2006l) . who 
used the technique of spectroastrometry to show the Ha emis- 
sion is in fact due to the close companion to /? Cep. This com- 
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panion may well be a classical Be star and if this is the case the 
contradictory nature of the system to current Be paradigms is 
negated. Schnerr et al. (2006) used spectra in the R band about 
Ho-. Here we use spectroastrometry to probe the blue region of 
the companion's spectrum to test the hypothesis that the close 
companion is a classical Be star 

Spectroastrometry is a technique which utilises the spatial 
information present in a longslit spectrum. The information is 
contained in the spatial profile of the spectrum, specifically in the 
photocentre centroid and the spectral profile's width. Changes in 
the flux distribution as a function of wavelength are manifest 
by changes in the centroid and width of the spectrum. An un- 
resolved binary system with one star dominating the flux at an 
emission line is revealed by a centroidal displacement towards 
the dominating star over this line. Conversely, if one star of the 
system has a strong absorption line in its spectrum the centroid 
of the spectral profile will shift to the other star over this line. 
Such signatures can be detected with high precision, of the or- 
der of 1 mas or less jOudmaiier et al. I l2008h . Therefore this is a 
powerful technique wi th which to detect and study close binary 
systems, as shown bv Baines et al.l (I2006.) . 

The spectroastrometric signature of an urrresolved binary 
system contains information on the distribution of the flux ema- 
nating from the system. Thus spectroastrometry can be not only 
be used to detect binary systems, it can also deconvolve the ob- 
served spectrum into the individual spectra of its components. 
Here we use this technique to disentangle the spectra of the 
close p Cep binary components to investigate the properties of 
the close companion. 

This paper is structured as follows: in Sect. |2] the details of 
the observations are presented, in Sect. |3] we present our spec- 
troastrometric results, in Sect. |4] we discuss the spectra splitting 
methods and in Sect. |5] we present the results of splitting the 
spectra. A discussion and interpretation of the results follows in 
Sect.|6]and finally in Sect.|2]we summarise our findings. 

2. Observations and data reduction 

2.1. Observations 

The data presented here were obtained using the 4.2m William 
Herschel Telescope (WHT) with the Intermediate Dispersion 
Spectrograph and Imaging System (ISIS) spectrograph. The data 
were obtained on the 7th of October 2006. Spectra of the /3 Cep 
system in the B (4200-5000A) and R (6200-6900A) bands were 
taken simultaneously using the dichoric slide of ISIS. The slit 
width was to set to 5". The R1200 and B1200 gratings were 
used and the resolving power of the spectrograph was approxi- 
mately 3800 (measured from telluric lines). The Marconi2 and 
EEV12 CCDs were used on the red and blue arm respectively, 
each with a pixel size of 13.5;um. This resulted in angular pixel 
scales of 0.20"in the blue and 0.22"in the red region. As the 
average seeing during our observations was 1.27"the spectral 
profile was well s ampled, which is a requirement for accurate 
spectroastrometry (lBailevl[r9 98). 

The data were gathered as part of wider study of binary sys- 
tems. A wide slit (5") was used to ensure all the light from a 
given system entered the slit, despite the effect this had on the 
spectral resolution. Multiple spectra were taken at the follow- 
ing position angles (PA) on the sky: 0°, 42°, 68°, 90°, 180° and 

o o o I 1 

270. Data were taken at a PA of 42 and 68 as ISchnerr et alj 
( |2006|) suggested the binary was orientated at ^4-2 . Dispersion 
calibration arcs were made using CuNe and CuAr lamps. 



2.2. Data reduction 

Data reduction was conducted using the Image Reduction and 
Analysis Facility (IRApfl and routines written in Interactive 
Data Language (IDL). Flat field and bias frames were combined 
and the averaged flat field was then normalised. The raw data 
were then corrected using the averaged bias frame and the nor- 
malised average flat frame. Saturated exposures were discarded. 
The total intensity longslit spectra were then extracted from the 
corrected data in a standard fashion. Wavelength calibration was 
conducted using the arc spectra taken after the science observa- 

o 

tions at a position angle of 270 . 

Spectroastrometry was performed by fitting Gaussian pro- 
files to the spatial profile of the longslit spectra at each disper- 
sion pixel. Spurious fits (for example due to cosmic rays) were 
identified and discarded, allowing the routine to fit the spectral 
profile. This resulted in a positional spectrum - the centroid of 
the Gaussian as a function of wavelength - and a FWHM (Full- 
Width-at-Half-Maximum) spectrum - the FWHM as a function 
of wavelength. The continuum position exhibited a general trend 
across the CCD chip (of the order of 10 pixels). This was re- 
moved by fitting a low order polynomial to the continuum re- 
gions of the spectrum. This set the continuum position of the 
centroid to zero. Spot checks were used to ensure line effects 
were not fit by the function. 

A correction for slight changes in the dispersion was deter- 
mined by cross correlating individual intensity spectra, and was 
then applied to the associated intensity, positional and FWHM 
spectra. This was to ensure slight changes in wavelength (due 
to flexure of the spectrograph) did not introduce spurious sig- 
natures when spectra obtained at differing position angles were 
combined. All intensity, positional and FWHM spectra at a given 
position angle were then combined to make an average spectrum 
for each position angle. 

The average positional spectra for anti-parallel position an- 
gles were combined to form the average North-South (NS) 

o o 

and East- West (EW) positional spectra, i.e.: (0 - 180 )/2 and 

o o 

(90 - 270 )/2. This procedure eliminates instrumental artifacts 

o 

as real signatures rotate by 180 when viewed at the anti-parallel 
position angle while artifacts remain at a constant orientation. 
However, in some of the data a signature was noted at only one 
position angle. Subtraction of the two anti -parallel spectra would 
not remove this effect but clearly it is also an artifact. Thus all av- 
erage position spectra were assessed visually to exclude features 
only present at a single position angle. Signatures in data taken 
at position angles with no anti-parallel counterpart were judged 
real if they exhibited qualitative similarity to signatures at a sim- 
ilar position angle that were judged real. The FWHM spectra at 
opposite position angles were visually inspected to search for ar- 
tifacts. Features present in data taken at only one position angle 
were classified as an artifact and discarded. FWHM spectra at 
anti-parallel position angles were combined to make an average 
NS or EW FWHM spectrum. 



' IRAF: written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO) in Tuscon 
Arizona ( h ttp://iraf.noao.edu/ 1) 
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3. Results: The spectroastrometric signatures 
observed 

3.1. The spectrum of/3 Cep 

In Fig. [T]we present the average spectra observed in both the B 
and R range. Ha can be seen in absorption with a small double 
peaked emission profile, with the blue peak just rising above the 
continuum level. The hydrogen Balmer line profiles are also pre- 
sented in the top panels of Fig. |2l No emission is noticeable in 
the Hy6 and Hy lines. Besides Ha other prominent lines in the red 
region are Hei at 6678A and a Diffuse Interstellar Band (DIB) 
at ^6280A. In the blue region many more absorption Unes are 
present, most prominent of which are the lines due to H i and 
He I. The numerous weaker and narrow lines present are primar- 
ily due to 11, Si ii. Si iii and other ionised metals such as C ii and 
Nn. 

3.2. The spectroastrometric results 

The spectroastrometric signatures of the close j8 Cep binary sys- 
tem over the three principle hydrogen Balmer lines are presented 
in Fig.|2] For each line we present: the intensity profile and the 
centroidal signature and change in the FWHM of the longslit 
spectra over these lines in the NS and EW directions. The photo- 
centre of the spectral profile shifts to the North and East and the 
FWHM increases over these lines. This means we clearly detect 
the close binary system - something not trivial in seeing condi- 
tions of greater than 1 " . The signature is most prominent over the 
Ha line while small and consistent features can be noted across 
Hj3 and Hy. The spectroastrometric signature observed over the 
emission component of the Ha profile implies that the source of 
the emission lies to the North-East (NE) and the emission profile 
is intrinsicall y broader than the ab sorption profile. This confirms 
the result of Schnerr et al.l (12006) . i.e. the Ha emission is ema- 
nating from the companion in the NE and not the primary. 

Across the H/3 and Hy lines the positional excursions occur 
in the same direction as over Ha. As these signatures take place 
across an absorption line this indicates that the South- West (S W) 
component of the system 'dominates' the absorption profile. In 
addition narrow positional excursions to the NE and FWHM in- 
creases can be seen over the lines seen alongside Hy (primarily 
On lines). This indicates that these lines are associated with the 
SW component of the system, the primary. These data highlight 
the exquisite sensitivity of spectroastrometry to changes in the 
flux distribution. Changes in the photocentre of the order of 5 
mas or less are detected. The noise in the positional spectra is 
typically of the order of 1 mas. 

The 'XY plots' (NS against EW excursions) of the /3 Cep 
system are presented in Fig. [3] The XY plots trace a straight line 
to the North and East - the direction in which the source of the 
Ha emission lies. The position angle of the binary is determined 
from a simple least-squares fit to the data, and is in general con- 
sistent across the different lines (Ha: 51.3 + 0.7°, Hfi: 66.5 ±4.6° 
& Hy: 58.7 ± 6.7°). However, the H/3 position angle is not con- 
sistent with the Ha data. The Hy6 and Hy lines were associated 
with much weaker positional excursions and as such these po- 
sition angles are more uncertain than the position angle derived 
from the Ha data. 

4. Splitting the spectra 

Two approaches were used to split entangled/c omposit e binar y 
spectra. The first approach was pioneered by iBailevI (Il998h . 



This method utilises the fact that the changes in the flux ra- 
tio of a binary system lead to positional displacements of the 
photocentre. The movement of the photocentre across a given 
wavelength is proportional to the component separation and the 
component flux ratio at the particular wavelength. This centroid 
movement takes place about the continuum position, which is 
determined by the continuum brightness ratio of the two compo- 
nents. Therefore if the separation and the continuum flux ratio 
of the two components are known the intensity spectra and po- 
sitional spectra observed can used to disentangle the individual 
fluxes of the two components. 

The second approach, that of iPorter et alj (|2004|) . does not 
require any prior knowledge of the binary system. This method 
not only deconvolves spectra, it also estimates the separation of 
the two binary components. This method is based upon numer- 
ical simulat i ons of point source binary systems performed by 
iPorter et al] ( |2004 . iPorter etalj (|2004 determined the depen- 
dence of the spectroastrometric signature of a given binary on the 
system's properties. Using the relationships established by these 
simulations one can use the three spectroastrometric observables 
(the centroid, total flux and width at a given A), with knowledge 
of the seeing, to recover the individual fluxes of the binary com- 
ponents. For a full description of the method see IPorter et aU 
(2004). 

We have used both methods for t wo reasons. Firstly there 
are certain limitations to the method of 'P orter et al.l (12004 '). We 
found accurate knowledge of the seeing was essential, something 
not trivial in this case. To test the reliability of the method of 
IPorter et al] 12004) we applied this method on simulated data. 
It was found that with accurate knowledge of the seeing this 
method consistently returned the separation of the two point 
sources and their respective fluxes, to within an error of a; 10% or 
better However, when adopting the minimum FWHM as an esti- 
mate of the seeing (which is actually an upper limit) the method 
can under estimate the separation. While we expect this effect 
to be small we treat our estimate of cr^ with caution and use the 
method of Bailev ( 1998) as a consistency check. 

Secondly, we do not just use the method of IBailevI (Il998h 
as although we can find literature values for the difference in 
brightness of the two components and their separation, none are 
certain. The difference in magnitude reported ranges from 1.8 
to 5.0 tBalegaetalJ 120021: iGezariet all 1 19721) while the sepa- 
ration changes slowly due to the motion of the binary system . 
The orbit has been determined by P igulski & Boratvnl ( Il992h . 
H owever this orbit is n ot consistent with our results and those 
of lSchnerr et al.l (l2006l) ( see Sect.l6.2li. There f ore, w e do not use 
the orbital parameters of iPigulski & Boratvnl d 19921) to estimate 
the position of the secondary. Thus we used each method as a 
consistency check upon the other. 

4.1. The method of Bailey (1998) 

The method of iBaileyl (Il998l) uses two spectroastrometric ob- 
servables (centroidal excursions and the total flux observed) and 
two inputs (binary projected separation and difference in bright- 
ness). Three values for the difference in brightness between the 
binary components were taken from the literature: 1.82 magni- 
tudes at A ^ 8100A, a rough average value of 3.4 magnitudes 
around 5500A and 5 magnitudes at 5000A dBalega et al.ll2002l: 
iHartkopf et aIll200ll:lGezari et al.lll972 '). Distances in the litera- 
ture between the primary and the secondary range from kQ.25" 
to 0.04". The most recent observation of the system was in 1998, 
and thus the position of the secondary in 2006 is uncertain. 
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Fig. 1. The average spectra of yS Cep (total light of the binary system) in the blue (upper panel) and the red (up to the telluric lines - 
lower panel) spectral range. 
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Fig. 2. The spectroastrometric signature of the close fi Cep binary system across the Ha, and Hy lines (from left to right). Across 
each line we present the averaged intensity profile normalised to the continuum, the position of the photocentre with respect to the 
continuum position in the NS and EW direction and the FWHM of the spectral spatial profile, with respect to the continuum value, 
in the same directions. In the positional spectra North and East are presented as being in the positive direction (i.e. to the top of the 
page). 



We used a range of distances (0.05"to 0.25") with each of 
the magnitude differences listed above, and evaluated the results 
based on the following criteria. Any input parameters leading 
to negative flux in the secondary spectra were immediately dis- 
regarded. Also situations where absorption lines in the primary 
were mirrored by an emission line in the secondary spectra were 
considered unlikely. 

In the red region a brightness difference of 5 magnitudes was 
discarded as it resulted in secondary spectra with negative flux. 
A difference in brightness of 3.4 magnitudes was discounted 
as this led to secondary spectra with Hei 6678A in emission, 
an uncommon occurrence in field stars. Using A/? of 1.8 with 
a range of separations did not result in a significant constraint 
upon the separation of the two components. A separation of 
0.05"was discarded as this resulted in Hei emission in the sec- 
ondary spectrum. Separations of 0.1"to 0.15"resulted in a sec- 



ondary spectrum devoid of any He i 6678 A feature while sepa- 
rations of 0.20"and 0.25"led to a secondary spectra with a He i 
6678A absorption feature. If the secondary is a mid B type star 
the presence of He i 6678A is to be expected, but the absence 
of a He I feature is consistent with a late type B star. Both spec- 
tral types are possible and thus the binary separation is not con- 
strained beyond the range 0.05"and 0.25". 

In Fig. |4] we present the average red spect ra of the clos e 
P Cep binary system, split using the method of lBailevI ( Il998l) . 
a brightness difference of 1.8 magnitudes and separations of 
0.05", 0.15" and 0.25". It is evident that the Ho- emission is 
associated with the secondary. The double peaked Ha emission 
profile is typical of rapidly rotating classical Be stars. 

When splitting the blue spectra differences in brightness be- 
tween the two binary components of 5.0 and 3.4 were discarded 
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Fig. 3. XY plots of the spectroastrometric signature of the close p Cep binary system over the Ha, Hj6 and Hy lines. The solid line 
is a fit to the Ha data, the dotted line is a fit to the Hy6 data and the dashed line a fit to the Hy data. North is to the top of the page 
and East is to the left. 




1.0 
0.8 
0.6 




Primary ^'''^\:/^ 
■ . 1 . . r . 1 . . . . 1 . . . . 1 . . . 


1.6 

1.2 


. . 1 , , , , , . 
" Secondary 


1 , , . . 

/" \ _ _ 




0.8 


■ . 1 . . . . 1 . 







6540 6555 6570 6585 
Wavelength (A) 



6670 6675 6680 6685 
Wavelength (A) 



1.2 
J 0.8 
S 0.4 

0) 3 
^ 2 



— Primary 

H 1 1 1 H 



Ha 

\ 

H 1 1 



He I 



Secondary 



_i_ 



+ 



TV 

u — 



0) 
E- 



J I I I I I I _i I I I I I ■ 



6300 6400 6500 6600 

Wavelength (A) 



6700 



6800 



Fig. 4. The red spectrum of the close /3 Cep binary system, split using the method iBaileyl (Il998l) . AR of 1.8 and a range of binary 
separations: 0.05 "(long-dashed line), 0. 15"(short-dashed line) and 0.25 "(solid line). In the top left the wavelength range is reduced 
and centred on Ha while in the top right attention is focused upon the He i 6678. 15A profile. 



due to similar arguments as presented above. Thus it appears the 
brightness difference between the two binary components in the 
B band is approximately 1.8 magnitudes, as in the R band. In 
Fig.|5]we present the spectra obtained with this value of AB and 
a separation of 0.20". The many narrow absorption lines present 
in the composite spectrum (On, Sim etc.) are clearly associated 
with the primary spectrum. Also evident is that the He i absorp- 
tion lines are weaker in the secondary spectrum than in the pri- 
mary spectrum. The secondary spectra obtained using smaller 
separations were judged to be of dubious validity as they exhib- 
ited many emission features coincident with absorption features 
in the primary spectrum. 

4.2. The method of Porter et al (2004) 

The model of iPorter etal] ( 12004 ') is based on a constant con- 
volving function, yet a change in focus along the chip length 
led to a change in the width of the flux distribution. To negate 
this difficulty the cr distribution was normalised via a polyno- 



mial fit to remove changes not due to the binary system. The 
mean seeing was estimated from the minimum of the now flat- 
tened FWHM spectrum. However, this estimate of the seeing 
is an upper estimate. Simulating a binary system with a sepa- 
ration of 0.20"and a difference in brightness of 1.8 magnitudes 
we found the minimum FWHM may over estimate the seeing 
by 0.02 pixels. Thus we subtracted this from our estimate of the 
seeing using the FWHM minimum. 

In figure |6] we present th e red spectrum of the j3 Cep sys- 
tem split using the method of iPorter et al.l (|2004|) . The spectrum 
exhibits a close similarity to the spectrum obtained using the 
method of Bailey ( 1998) and similar parameters as those output 
by the method of IPorter et all (120041) . 

Applying the method of IPorter et al.l (l2004t) to the blue spec- 
tra did not result in separated spectra, as the model did not con- 
verge upon the most likely separation. This could be due to a 
lack of prominent features in the position spectra. The precision 
of the blue position spectra was no worse than that in the red 
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Fig. 5. The blue spectrum of the close p Cep binary system, split using the method iBailevI (1 1998 *). A B = 1.8 and d - 0.20". In the 
top left the wavelength range is reduced and centred on Hy while in the top right attention is focused upon the Hy6 profile. 



Primary 



0) 




6540 6555 6570 6585 
Wavelength (A) 



6670 6675 6680 6685 
Wavelength (A) 




6300 6400 6500 6600 

Wavelength (A) 



6700 



6800 



Fig. 6. The red spectrum of the close /? Cep binary system, split using the method [Porter et all (|2004|) . In the top left the wavelength 
range is reduced and centred on Ha while in the top right attention is focused upon the He 1 6678. 15A profile. The determined values 
for the separation and brightness difference were d-Q.2\ + 0.07" and A/?=1.5 + 0.17 



spectra (of the order of 1 mas). However, the positional offsets 
in the blue region are much smaller than in those in the red re- 
gion as the flux difference in the blue region is less pronounced 
than that over the Ha line. 



5. Results: the separated spectra 

5.1. The spectral types of the binary components 

The spectral type of /? Cep has been estimated to be both B 1 and 
B2 while the luminosity class of the star has been report ed to be 
III, ly and V (Morgan et al. 1943, 1955; Lesh 1968; More l et alj 
|2006|) . These spectral classifications refer to the total light of the 
system. As we have disentangled the spectra of the binary com- 



ponents we can now spectrally type each component separately. 
The spectrum of the binary system has also been designated 'e 
& v' where e refers to the emission nature and v the variability 
of the spectrum. We assume the that the primary is responsible 
for the variability of the spectrum due to its pulsating nature. To 
determine the spectral type of the binary components w e use the 
spectra that were separated with the method of Baile^ (Il998l) . a 
brightness difference of 1 . 8 magnitudes and a separation of 0.2" . 

The presence of He i absorption in the red spectrum of the 
primary indicates a B type star. The equivalent widths of the Ha 
and He i absorption features suggest that the primary is an early 
B (1-3) type star, with a luminosity class of IV or III. In the 
blue region the primary spectrum exhibits prominent H i absorp- 
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tion lines (H/3 and Hy) alongside Hei absorption lines, which 
indicates a B type star. The relatively strong He i lines and com- 
paratively weak H i lines indicate that the primary is an early B 
type star (i.e. BO to B2). Additional features in the primary spec- 
trum include absorption lines due to He i, O ii, Si iii, C ii, C iii and 
N II, lines which are also suggestive of an early B type star No 
He II lines are present, implying the star has a spectral type later 
than 09. The ratio of the He 1 447 1 A/Mg ii 448 1 A lines indicate 
the star is an early B type star, most probably B2. To determine 
the luminosity class of the star the ratios of the On 4415 A & 
4417A to He 1 4387A and O ii 4348A to Hy lines were used. The 
strength of the On lines indicates that the primary's luminosity 
class is most probably III. 

To summarise, the spectral type of the primary was deter- 
mined to be B2III. The uncertainty associated with this conclu- 
sion is approximately one spectral subtype. 

The most prominent feature in the red spectrum of the sec- 
ondary is that of the Ha emission. The Ha emission extends 
approximately +400 km s~' f rom th e rest wavelength of Ha, as 
also found by ISchnerr et alj (l2006h . A Gaussian fit to the Ha 
emission profile results in a value for the FWHM of the profile 
of a! 500km s '. This width is typical of a rapidly rotating Be 
star. A He i 6678 A absorption feature is present, although weak, 
with an equivalent width of * 0.25 A. This places an upper limit 
on the spectral type of about B4/B5. 

The secondary spectrum in the blue does not display the 
many narrow absorption features due to O ii and Si iii the primary 
does. Thus there is no reason to suspect it is not a dwarf star. The 
spectrum of the secondary in the blue region exhibits absorption 
lines due to H i {Hfi and Hy) and He i (i.e. 4388, 447 1 A). This is 
indicative of a B type star. The He i lines are weaker than than 
those in the spectrum of the primary (equivalent widths are ap- 
proximately 15% less) which suggests the secondary is of a later 
spectral type than the primary. This is to be expected as it is less 
bright. In contrast, the H i lines in the secondary are of a similar 
strength to those in the primary, indicating an early B type star 
However, as the secondary shows emission in Ha it is possible 
there is some emission in Hfi and Hy which 'fills in' the absorp- 
tion profile of the secondary. Concentrating on the strength of 
the He I lines, and the ratio of Hei/Mgii (447 1/4481 A) fines, we 
conclude the most likely spectral type of the secondary is B5. 

Therefore the spectral type of the secondary is determined to 
be B5 Ve. The uncertainty in this spectral typing is again approx- 
imately one spectral subtype. 

In determining the spectral type of the secondary we assume 
that all the flux observed emanates directly from the star in ques- 
tion and the spectral features observed are purely photospheric in 
their origin. To asses the possible flux c ontribution from circum - 
stellar material we use the NIR study of lDoughertv et al] (Il99ll) . 
In most cases mid type Be stars were found to have a small V- 
J excess, the ave rage (V-J) excess of a B5Ve star according to 
[Dougherty et al.l ^1991) is only ~6% of the stellar continuum. 
In addition the continuum excess falls off very sharply with de- 
creasing wavelengths. Therefore we expect the continuum ex- 
cess due to any circumstellar material to be negligible at optical 
wavelengths. Regarding emission lines there may well be unre- 
solved H I lines in the secondary spectrum therefore we do not 
rely on the strength of the H i lines to determine the spectral type 
of the secondary. There may also be some emission component 
in the He i lines. Howeve r, in a s tudy of more than 100 Be stars 
over a period of 10 vears lChauville et al.l (1200 ID did not observe 
any Be star with net He 1 4471 A emission. Therefore, in light of 



the above comments, we assume that the He i lines observed are 
photospheric in origin. 

5.2. The mass ratio of the close binary system 

Taking the spectral type of the primary to be B2III and that 
of the secondary to be B5 V and using th e tabulated values of 
iLandolt et alj ( IT98Z) and Harmanedd 19881) we obtain values for 
the masses of the components of the system of Mpri=12.6 + 
3.2M0 and Msec=4.4 + O.7M0 respectively. Uncertainties were 
estimated by allowing an uncertainty in the spectral types of 
one subtype. The mass ratio of the system is determined to be 
0.35 + 0.15. 

5.3. The V sin / of the individual binary components 

Using the separated spectra the v sin / values of each component 
of the binary system can be estimated. A rotational profile was 
constructed for a variety of vsin/ values: from 1 to 550 kms"' 
in steps of 1 kms The rotation profile was then convolved 
with synthetic spectra constructed using ATLAS9 and S YNTHE 
(lKurucdl993lR Following the convolution of the rotational pro- 
file and the synthetic spectra the resultant spectra were broad- 
ened by convolution with a Gaussian function to match the spec- 
tral resolution. The rotationally broadened profile of the Hei 
447 1 A line was then compared with the observed profile using 
a simple test. 

For the primary synthetic spectra were constructed using val- 
ues of Teff of 23000K, 24000K and 25000K and lo g(g) of 3.8 
These values were based on the values determined bv lCatanzarol 
(2008) and the spectral type determined previously. To simulate 
the spectrum of the secondary synthetic spectra were constructed 
using values of T^s of 14000K, 15000K, 16000K and 17000K 
and log(^) of 4.1. These v alues are based o n a B5 type dwarf 
star on the Main Sequence dHarmaned 19881) . The range of tem- 
peratures used reflects the uncertainty in the spectral typing. A 
micro-turbulence of 2 kms"' was used in the synthetic spectra 
generation, and when constructing the rotational profiles a limb- 
darkening coefficient of 0.6 was assumed. 

The best fit v sin / values obtained are dependent upon the 
value of Teff used to generate the synthetic spectrum. This is be- 
cause the depth of the Hei line increases towards earlier types. 
With the previous spectral typing an error in T^ff of » 1 OOOK pos- 
sible. However, the input T^ff does not effect the width of the 
convolved profile. Thus the most likely value of Teff is the one 
for which the rotationally broadened synthetic spectrum matches 
not only the height but also the width of the observed profile 
at a given value of vsin/ (i.e. that associated with the smallest 
X^). This effect is illustrated in Fig.|2]by plotting the best fit ro- 
tationally broadened profiles at 15000K and 17000K over the 
secondary profile. The best fit at 15000K fails to fit the width 
and depth of the observed profile, which was reflected in the 
large value for the minimum at this temperature. The best fit 

V sin / values were lower for lower values of Teff and ranged from 

V sin /=132 km s"' at Teff of 14000K to v sin /=230 km 5"' at Teff 
of 17000K. The lowest values were obtained with a T^fi of 
24000K in the case of the primary and 17000K in the case of the 
secondary. We determine that the primary has a v sin / of 29;^29 

- We used the GNU Linux port of ATLAS9 and SYNTHE developed 
bv lSbordoneetal.1 ([2004). For the initia l atmo spheric models we used 
the grid of models bV iCastelli & Kuruczl ( l2003h . 
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Fig. 7. A comparison between the spectroastrometrically sep- 
arated He I 4471 A profiles and the best fit rotationally broad- 
ened synthetic profiles. The solid line represents the Hei line of 
the primary and secondary spectra respectively, the short-long 
dashed line the best fit primary synthetic spectrum (Teff 24000K, 
\ogig) 3.8, vsin/ 29 kms"'), the short-dashed line the best fit 
secondary synthetic spectrum at a T^^ of 15000 K {\og{g) 4.1, 
vsin / 163 kms ') and the long dashed line is the best fit sec- 
ondary synthetic spectrum at a T^g of 17000 K (log(g) 4. 1, v sin / 
230 kms-'). 



kms^ and the secondary has a vsin / of 230 + 45 kms . The 
uncertainties in these values were determined from the change 
in V sin / of the synthetic spectra which led to an increase in the 
of the fit of 1 . We note that the secondary profile was best fit 
using a B4 star profile {T^-s - 17000K), rather than a B5 profile. 
This is consistent with an uncertainty of 1 subtype in the spectral 
type of the secondary. 

6. Discussion 

6.1. The V sin / of the individual binary components 

The V sin / of the primary is consisten t with literature values of 
around 25 kms ' dTelting et al.lll997h . although it is imprecise. 
This imprecision is due to the use of a relatively wide slit result- 
ing in an instrumental profile with a FWHM of ^10 km s '. The 
vsin / value for the secondary is similarly imprecise. However, 
it is clear that the secondary does rotate substantially faster than 
the primary, in agreement with the estimate of si230 kms"' of 
[Catanzaro (2008) . 

Given the spectral type of t he secondary it's critic al velocity 
is approximately a;430 km s~' dTownsend et al.ll2004 V It is sug- 
gested that the o rbit of the yS Cep system is seen almost edge on 
JPigulski & BoratvnII 19921 see also Sect. 6.2). If the stars rotate 
in the plane of the orbit, which is not necessarily the case, then 
the above v sin / values are the intrinsic rotation velocities of the 
stars. In this case the rotational velocity of the secondary is 45- 
65% of its critical velocity. This rotation may not be consistent 
with the hypothesis that the secondary is a classical Be star. 

Whether or not all Be stars rotate at their break-up veloc- 
ity is currently a matter o r some deba te, with argum ents both 
for (iTownsend et al.ll2004l) and against ( ICranmeij|2005h this sce- 
nario. However, it is generally accepted that all Be stars rotate 



at a substantial fraction of their break-up velocity. To prove the 
secondary is rotating at/near its critical velocity we need an inde- 
pendent determination of its inclination, which is far from triv- 
ial. In addition even if / was constrained, this v sin / measurement 
would not be conclusive. It has been shown the method used here 
to determine v sin / re turns a lower Uniit as it does not account 
for gravity darkening (iTownsend et al.ll2004 ). However, the sec- 
ondary is shown to be rotating relatively fast, even at the lower 
limit of / = 90°, and it may be rotating at a substantial fraction 
of its bre ak-up v elocity. This is essentially consistent with the 
finding of lPortej ([1996) who demonstrated Be stars intrinsically 
rotate at approximately 70% of their break-up velocity (using a 
similar technique to estimate v sin /). Therefore it is indeed hkely 
the secondary is a classical Be star. 

6.2. The orbit of/3 Cephei 

The orbit of the close companion to /3 Cep has been deter- 
mined bv lPigulski & Bo ratvn ( 1 9921) . Howe ver, both the results 
presented here and by Schnerr et al.l (12006 ) suggest a revision 
of t he orbital parameters i s required. According to the orbit of 
Pig ulski & Boratviil d 19921) the companion to /3 Cep should have 
been t o the S W of the pri mary in 2006. The data presented here, 
and bv lSchnerr et al.l (|2006) are not consistent with this predic- 
tion. As the spectroastrometric data proves the Ha emission is 
emanating from the close companion, the data place the com- 
panion to the NE of the primary in 2006. 

Here we investigate whether a slight change of orbital pa- 
rameters can result in an orbit consistent with observations. It 
was found the data allow a straight line fit, which implies that 

o 

the system is viewed at a very high inclination, i.e. a!90 . From 
a least-squared fit to the data we obtain a value for the position 

o 

angle of the line of nodes to be 228.6 + 1 . 4 . We took the date 
of 1914.6 as our reference periastron, as did Pigulski & Boratyn] 
(.1992). For the period of the orb it we considered the suggestion 
of iHadrava & Harmaned d 19961) that the system was approach- 
ing periastron in 1996. This period of 81.4 years differs from 
the value used by iPigulski & Boratvnl d 19921) but is within 3cr 
of their value. The semi-major axis of the orbit was estimated 
via Kepler's third law, the above period and the total mass of 
the system. With the data at our disposal constraining e and o) 
is difficult. We set e at 0.6 and investigated what value of a> was 
required to fit the observational data, given the parameters deter- 

o 

mined above. We found a value of a> of 20.0 fit the data well 
(see Fig. [8]). 

The final orbital parameters used were: Porh - 81.4yr, 
e = 0.60, To = 1914.6, Q = 228.6°, i = 90.0°, a = 0.25" 

o 

and u) - 20 . Most parame ters are within 3cr of the values 
of tPigulski & Boratvnlll992l) . The orbit is consistent with pre- 
vious speckle interferometric observations of the system and 
our spectroastrometrically determined position of the secondary. 
Admittedly, systems are rarely observed at an inclination of ex- 

o 

actly 90 . Changing the inclination by up to a few degrees does 
little to qualitatively change the picture. Provided Q. is revised 
in light of any inclination changes, a consistent fit to the data 
is achieved when / is changed by a few degrees (see Fig[8]l. We 
stress such an orbit is only illustrative, and more observations are 
needed to fully constrain all the orbital parameters. However, we 
successfully demonstrate a slight revision of the orbital parame- 
ters of Pigulski & Bora tvn ( 1992) is all that is required to obtain 
an orbit consistent with observations. 
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Fig. 8. A possible relative orbit of the companion of the 
close yS Ceph system about the primary. Speckle interferomet- 
ric data have b een taken from Pigulski & Boratvn (1992) and 
iHartkopfetaf] (120011) . and are marked by crosses. The open 
square point represents the position of the companion from our 
spectroastrometric data, the line indicates the orbital path. We 
compare specific orbital predictions and data points at three 
epochs, 1971.5, 1998.8 and 2006.8. We denote these data points 
by open triangles (a square in the case of 2006.8) and orbital 
predictions by open circles. The primary is situated at 0,0 and 
the filled circle marks the position of periastron. The dashed 

o 

line shows the orbital path when / is changed by -2 and Q is 

o 

changed by -3 . The average uncertainty in the data was ap- 
proximately 0.015" in each direction. 



7. Conclusions 

In this paper we have studied the close binary system of p Cep 
utilising the novel approach of splitting the binary spectra us- 
ing spectroastrometry. Disentangling the composite binary spec- 
trum allows us to determine key properties of each component. 
We find that the Ha e mission of the systern is due to the close 
secondary, as shown bv lSchnerr et al.l ( l2006l) . Splitting the con- 
volved spectrum and assessing the spectral type of each compo- 
nent we find that the companion is a dwarf star with a spectral 
type of B5. We also find it may rotate at a substantial fraction 
of its critical velocity, with a lower limit of v sin / - 230 km s"' 
corresponding to V^IVan of 53%. 

The secondary's estimated mass and vsin/ fall within the 
range of typical values for classical Be stars. Thus we consider 
it highly likely the secondary is indeed such a star. Therefor e we 
have not only confirmed the result of ISchnerr et al.l (l2006h but 
we validate their suggestion that the secondary could be a clas- 
sical Be star. In which case the Ha emission is thought to be due 
to gaseous equatorial material ejected from the star by a com- 
bination of rapi d rotation and some othe r phenomena, e.g. non- 
radial pulsation ( Porter & Riviniusl2003h . As the Hot emission is 
shown to originate from a star where the standard Be paradigm 
applies the [5 Cep system does not pose any contradictions to the 
current understanding of either/? Cep stars or classical Be stars. 
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